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Abstract 

We consider atomic chains with nonlocal particle interactions and prove the existence of near- 
sonic solitary waves. Both our result and the general proof strategy are reminiscent of the seminal 
paper by Friesecke and Pego on the KdV limit of chains with nearest neighbor interactions but 
differ in the following two aspects: First, we allow for a wider class of atomic systems and must 
hence replace the distance profile by the velocity profile. Second, in the asymptotic analysis we 
avoid a detailed Fourier pole characterization of the nonlocal integral operators and employ the 
contraction mapping principle to solve the final fixed point problem. 

Keywords: asymptotic analysis, KdV limit of lattice waves, 

Hamiltonian lattices with nonlocal coupling 

MSC (2010): 37K60, 37K40, 74H10 

Contents 


1 Introduction 1 

1.1 Setting of the problem. 2 

1.2 Overview on the main result and the proof strategy. 3 

2 Preliminaries and linear integral operators 5 

2.1 Reformulation in terms of convolution operators . 5 

2.2 Asymptotic analysis for the convolution operators Ar] . 6 

2.3 Asymptotic properties of the auxiliary operator . 8 

3 Proof of the main result 10 

3.1 The leading order problem and the KdV wave. 11 

3.2 The linearized traveling wave equation for £ > 0. 13 

3.3 Nonlinear fixed point argument. 16 

List of symbols 19 

Acknowledgements 19 

Bibliography 19 


1 Introduction 

Since the pioneering paper [ZK65], the so-called KdV limit of atomic chains with nearest neighbor 
interactions - often called Fermi-Pasta-Ulam or FPU-type chains - has attracted a lot of interest 
in both the physics and the mathematics community, see [FML14] for a recent overview. The key 
observation is that in the limiting case of long-wave-length data with small amplitudes the dynamics 
of the nonlinear lattice system is governed by the Korteweg-de Vries (KdV) equation, which is 
a completely integrable PDF and hence well understood. For rigorous results concerning initial 
value problems we refer to [SWOO] and to [CBCPS12, GMWZ14] for similar result in chains with 
periodically varying masses. 
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Of particular interest are the existence of KdV-like solitary waves and their stability with respect 
to the FPU dynamics. Both problems have been investigated by Friesecke and Pego in the seminal 
four-paper series [FP99, FP02, FP04a, FP04b], see also [HW13] for simplifications in the stability 
proof and [FML14] concerning the existence of periodic KdV-type waves. The more general cases 
of two or finitely many solitary waves have been studied in [HW08, HW09] and [Mizll, Mizl3], 
respectively. In this paper we generalize the existence result from [FP99] and prove that chains 
with interactions between further than nearest-neighbors also admit KdV-type solitary waves. The 
corresponding stability problem is beyond the scope of this paper and left for future research. 


1.1 Setting of the problem 

We consider an infinite chain of identical particles which interact with up to M neighbors on both 
sides. Assuming unit mass, the equations of motion are therefore given by 

M 

ilj = ~ Uj) ~ ~ ; ( 1 ) 

m=l 

where Uj{t) denotes the position of particle j at time t. Moreover, the potential <hi describes the 
interactions between nearest-neighbors, 4*2 between the next-to-nearest-neighbors, and so on. 

A traveling wave is an exact solution to (1) which satisfies 

Uj{t) = r^j + v^t + sUeix), x:=ej-ecst, (2) 


where the parameters r* and u* denote the prescribed background strain and background velocity, 
respectively. Moreover, e > 0 is an additional scaling parameter which will be identified below and 
becomes small in the KdV limit. A direct computation reveals that the wave speed Cs as well as the 
rescaled wave profile Ug must solve the rescaled traveling wave equation 

M 

e^cl U” = ^ + me'^V+meUe) , (3) 

m=l 


where the discrete differential operators are defined by 


(V+™,V)(x) 


Y{x + me) — Y (x) 
me 


(V_^,K)(x) 


Y{x) — Y{x — me) 
me 


( 4 ) 


Note that u* does not appear in (3) due to the Galilean invariance of the problem and that the 
solution set is invariant under the addition of constants to Us- It is therefore natural to interpret (3) 
as an equation for the rescaled velocity profile Ws := U^; the corresponding distance or strain profile 
V+eC/g can then be computed by convoluting Ws with the rescaled indicator function of an interval, 
see formula (14) below. 

For M = 1 and fixed e > 0 there exist - depending on the properties of 4>i - many different types 
of traveling waves with periodic, homoclinic, heteroclinic, or even more complex shape of the profile 
Ws, see for instance [HerlO, HRIO, HMSZ13] and references therein. In the limit e —)■ 0, however, the 
most fundamental waves are periodic and solitary waves, for which Ws is either periodic or decays 
to 0 as X —)■ Too. 

In this paper we suppose r* = 0 - this condition can always be ensured by elementary transfor¬ 
mations - and split off both the linear and the quadratic terms from the force functions This 
reads 


<^'^{r) = amV + = O(r^) , m = 1...M (5) 

or, equivalently, ^m{r) = 5 amr^ + + 4'm(r) with 'l'm(r) = O(r^). In order to keep the 

presentation as simple as possible, we restrict our considerations to solitary waves - the case of 
periodic profiles can be studied along the same lines - and rely on the following standing assumption. 
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Figure 1: Sketch of the rescaled velocity profile Wg for e > 0 (black) and £ = 0 (gray) as function of the 
rescaled phase variable x. The grid with spacing £ > 0 describes the rescaled particle index ej while the dashed 
arrows indicate the height and the width of the pulse W^. The rescaled distance profile A^Wg has a similar 
shape. 


Assumption 1 (properties of the interaction potentials). For all m = the coefficients am 

and Pm are positive. Moreover, ^'m is continuously differentiable with 'l'm(O) = 0 and 

( 6 ) 

for some constants 7m and all r with |r| < 1. 

Note that the usual requirements for M = 1 are ai > 0 and Pi 0 but the case /3i < 0 can be 
traced back to the case /3i > 0 by a simple reflection argument with respect to the strain variable r. 
Below we discuss possible generalizations of Assumption 1 including cases in which the coefficients 
come with different signs. 

1.2 Overview on the main result and the proof strategy 

The overall strategy for proving the existence of KdV-type solitary waves in the lattice system (1) 
is similar to the approach in [FP99] but many aspects are different due to the nonlocal coupling. In 
particular, we base our analysis on the velocity profile 

m := U' (7) 

instead of the distance profile Veils, deviate in the justification of the key asymptotic estimates, 
and solve the final nonlinear corrector problem by the Banach fixed-point theorem. A more detailed 
comparison is given throughout the paper. 

As for the classical case M = 1, we prescribe a wave speed Cg that is slightly larger than the sound 
speed Co and construct profile functions that satisfy (3) and decay for x —)■ ±oo. More precisely, we 
set 


M 

cl := cl + £^ , cl := ^ aminnP > 0 , ( 8 ) 

m=l 

i.e., the small parameter e > 0 quantifies the supersonicity of the wave. Note that the subsonic case 
Cs < Co is also interesting but not related to solitary waves, see discussions at the end of §2 and the 
end of §3. 

The asymptotic analysis from §2 reveals that the limiting problem as e —)■ 0 is the nonlinear ODE 

IDo =diWo- d2W^ , (9) 

where the positive constants di and d 2 depend explicitly on the coefficient am and Pm, see formula 
(32) below. This equation admits a homoclinic solution, which is unique up to shifts (see §3.1) and 
provides via w{t, x) = ITo(x — t) a solitary wave to the KdV equation 


di dtw -b d2 dxw'^ -b d^w = 0 . 
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For e > 0 we start with the ansatz 

W, = Wo + eV, G LLn(K) (10) 

and derive in §3 the fixed point equation 

V,=Te[Ve] (11) 

for the corrector where the operator is introduced in (50). The definition of requires to 
invert a linear operator Ce, which is defined in (36) and admits a singular limit as e —)• 0. The linear 
leading order operator Cq stems from the linearization of (9) around the KDV wave Wq and can be 
inverted on the space Lgygjj(M) but not on L^(M) due to the shift invariance of the problem. The first 
technical issue in our perturbative existence proof is to show that this invertibility property persists 
for small e > 0, see Theorem 12. The second one is to guarantee that is contractive on some ball 
in Lgygjj(M), see Theorem 13. Our main findings are illustrated in Figure 1 and can be summarized 
as follows, see also Corollary 14. 

Main result. For any sufficiently small e > 0 there exists a unique even and nonnegative solution 
We to the rescaled traveling wave equation (3) with ( 8 ) such that 

\\W,-Wo\\2 + \\We-Wo\\^<Ce^ 

holds for some constant C independent of e, where Wq is the unique even solution to (9). 

The asymptotic analysis presented below can - for the price of more notational and technical effort 
“ be applied to a wider class of chains. Specifically, we expect that the following generalizations are 
feasible: 

1. We can allow for M = oo provided that the coefficients am, fm and 7m decay sufficiently fast 
with respect to m (say, exponentially). 

2. Some of the coefficients am and fm might even be negative. In this case, however, one has to 
ensure that the contributions from the negative coefficients are compensated by those from the 
positive ones. A first natural condition is 

M 

amm^ > 0 

m=l 

which ensures that uniform states are stable under small amplitude perturbations and that the 
sound speed cq from ( 8 ) is positive. A further minimal requirement is 

M M 

amm^ > 0 , fmm^ + 0 

m=l m=l 

because otherwise the leading order problem - see (9) and (32) below - degenerates and does 
not admit exponentially decaying homoclinic orbits. 

3. The non-quadratic contributions to the forces might be less regular in the sense of 

|^"(r)| 

for some constants 7 m and exponents 0 < Km < 1 - 

The paper is organized as follows; In §2 we introduce a family of convolution operators and refor¬ 
mulate (3) as an eigenvalue problem for W^- Afterwards we provide singular asymptotic expansions 
for a linear auxiliary operator B,., which is defined in ( 22 ) and plays a prominent role in our method. 
§3 is devoted to the proof of the existence theorem. We first study the leading order problem in 
§3.1 and show afterwards in §3.2 that the linear operator is invertible. In §3.3 we finally employ 
Banach’s contraction mapping principle to construct solutions to the nonlinear fixed problem (11) 
and conclude with a brief outlook. A list of all important symbols is given in the appendix. 
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2 Preliminaries and linear integral operators 


In this section we reformulate the nonlinear advance-delay-differential equation (3) as an integral 
equation and provide asymptotic estimates for the arising linear integral operators. 

2.1 Reformulation in terms of convolution operators 

For any rj > 0, we define the operator Ajj by 

1 fX+ri/2 

{Ar,Y){x) := - YiOdC (12) 

V Jx—ri/2 

and regard (3) as an equation for the rescaled velocity profile (7). Notice that Arj can be viewed as 
the convolution with the rescaled indicator function of the interval [—ri/2, +ri/2]. 

Lemma 2 (reformulation as nonlinear eigenvalue problem). Suppose that Ws = U'^ belongs to L^(M). 
Then, the nonlinear eigenvalue problem 


M 

e^clWe = ^ rnArns^'rni'IT^^'^-^rneWs) (13) 

m=l 

is equivalent to the traveling wave equation (3). 

Proof. The operators defined in (4) and (12) satisfy 

(y±meUe)ix) = {AmeUe)' {x ± |me) = {AmeWe) {x ± ^me) , (14) 

SO (3) follows from (13) after differentiation with respect to x and defining 14 as the primitive of 
Ws- In order to derive (13) from (3), we first notice that 114 £ L^(M) implies AmsWs £ W^’^(M) 
(cf. Corollary 4 below) and hence (AmeWs)(x) to 0 as x —)• ±oo. Afterwards we integrate (3) with 
respect to x and eliminate the constant of integration by means of the decay condition at infinity. □ 

In the case M = 1, we can derive from (13) the identity 

s^cjAsWe = Al^[{e^AeW,), 

which is the equation for the distance profile AeWg and has been studied in [FP99] (see equation 
(2.7) there for the function f = AeW^). For M > 1, however, we have to work with the velocity 
profile 114 since for a general function IF it is not possible to express AmeW for m > 1 in terms of 
AeW. 

We next summarize important properties of the convolution operators defined in (12). 

Lemma 3 (properties of Arj)- For each rj > 0, the integral operator Arj has the following properties: 

1. For any IF G L^(M), we have ArjW G n L°°(M) with 

\\A,W\\^ < r/-V2 liw^ll^ , \\A,W\\, < IIIFII 2 . (15) 

Moreover, A^W admits a weak derivative with ||(Flr;lT )^||2 < 2ri~^ IIITII 2 . 

2. For any IF G L°°(M), we have ||Fl,;lF||oo — Halloo- 

3. Arj respects the even-odd parity, the nonnegativity, and the unimodality of functions. The latter 
means monotonicity for both negative and positive arguments. 
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4- Ari diagonalizes in Fourier space and corresponds to the symbol function 


ar^{k) = sine {rjk/2) 


( 16 ) 


with sinc( 2 :) := sin( 2 :)/ 2 ;. 

5. Arj is self-adjoint in the \A-sense. 

Proof. All assertions follow immediately from the definition of Arj', see [HerlO] for the details. □ 

Corollary 4 (regularity of W G L^(]R) implies Ar^W G W^’^(M) C BC(M) and hence 

{AriW){x) —)• 0 as X ^ ±oo. 


2.2 Asymptotic analysis for the convolution operators 

The symbol function arj from (16) is analytic with respect to z = rik/2 and in view of 


1=0 




we readily verify 


A. 


- sinc(nk/2)e^^^ - - V 

^ - sme {r)k/2)e - Ij + 1)’ ~ 2^t{2j + 1)! ' 

The integral operator (12) therefore admits the formal expansion 


A = E 




221 (2j + 1)! ^ ^ ^ ’ 


1=0 


(17) 


which reveals that Ame should be regarded as a singular perturbation of the identity operator id. 
This singular nature complicates the analysis because the error terms in (17) can only be bounded 
in terms of higher derivatives. 

One key observation for dealing with the limit e —)■ 0 is - roughly speaking - that the resolvent- 
type operator 


^id+K 


id 


-1 


is well-defined and almost compact as long as k > 0. It thus exhibits nice regularizing properties 
which allows us to compensate bad terms stemming from the expansion (17). The same idea has 
been employed in [FP99] in the context of the distance profile AsW, showing that the Yosida-type 
operator 


^id -\-K 


id 


-1 


A 


2 

£ 


behaves nicely since the corresponding Fourier symbol 

e^a‘^{k) 

£2 + k(1 — a‘^{k)) 

is well-defined and bounded by (7/(1 -|- e^/c^), cf. [FP99, Corollary 3.4.]. Before we establish a related 
but weaker result in next subsection, we derive explicit error bounds for the singular expansion of 
Ame ■ 
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Figure 2: Left panel: Graph of the sine function z i—)■ sin(z)/z. Right panel Lower bound for 1 — sinc^ as 
used in the proof of Lemma 6 . 


Lemma 5 (small-parameter asymptotics of A-q ). There exists a constant C, which does not depend 
on rj, such that the estimates 


\\AqW -W\\^<Cif\\W" 


12 ’ 


\\AqW-W\\^<Cr,^\W'‘ 


and 


2 

AqW -W- ^W" 


< W'' 


12 ’ 


2 

A„W -W- ^W" 


24 


<Cp^ W'‘ 


hold for any sufficiently regular W. In particular, we have 


A„W > W strongly in V? 


(18) 


(19) 


( 20 ) 


for any W £ L^(M). 

Proof. -estimates : For any W £ W^’°°(M), the weak variant of Taylor’s expansion theorem implies 

(x - 0" 


P{x, 0 


< W" 


loo 24 

for almost all x, ^ G M, where 

P{x, C) := W{ff) - lF(x) - lT'(x)(x - ^) - \W''{x){x - if - lW'''{x){x - if . 
Integrating P{x, i) with respect to ^ G [x — ri/2, x -|- ry/2] we therefore get 


riAqW{x) - r]W{x) - f-,W"{x) 


24 


< 


/ Pn{x,i)di 
J x—rjA 

" / {x-ifdi = C\\W‘ 

Jx—n/2 


24 


w/ll 5 

v~i ' 


and ( 19)2 follows immediately. The derivation of ( 18)2 is similar. 

\A-estimates : Now let W £ W^’^(M) be arbitrary. By Parseval’s Theorem - and employing that 
|l — sine (z) — 2 :^/ 6 | < Cz^ holds for some constant C and all 2 : G M - we find 


2 

AqW -W -^W” 


W-Xjw 


1 — sine {r]k/2) — 


rj- 




24 


W{k) 


dk 


< Crf 


ffiW{k) dk = Crf W" 


12 ’ 


and this implies (19)i. The estimate (18)i can by proven analogously since we have |1 — sine {z)\ < 
z^ jd for all z G M. 
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Final argument: Let W G L^(M) be arbitrary but fixed. Since Arj is self-adjoint, see Lemma 3, 
and in view of (18) we readily demonstrate 


ArfW > W weakly in L^(M), (21) 

and this implies ||hL ||2 < liminfr,_ 5 .o ||^,jPL|| 2 . On the other hand, the estimate (15)2 ensures that 
limsup^_,.Q ||^,,VL ||2 < ||hL|| 2 - We therefore have ||VL ||2 = lim^_^o ll“^r?W ||2 and combining this with 
the weak convergence (21) we arrive at (20) since L^(M) is a Hilbert space. □ 


2.3 Asymptotic properties of the auxiliary operator 

As already outlined above, we introduce for any given e > 0 the operator 

M . , ^2 

Be ■=A+Y, , ( 22 ) 

m=l ^ 


which appears in (13) if we collect all linear terms on the left hand side, insert the wave-speed scaling 
(8), and divide the equation by We further dehne the operator 


Ho := id 


Z^m=l 




12 


dl 


(23) 


which can - thanks to Lemma 5 - be regarded as the formal limit of Be as, e ^ 0. In Fourier space, 
these operators correspond to the symbol functions 


M 


he{k) = l+Y. 


CXm^ 


2 1 — sinc^ {mke/2) 


bo{k) = 1 + 


m=l 


\-^M 

/^m=l 

12 


,m 


(24) 


which are illustrated in Figure 3 and satisfy 


be{k) ^ bo{k) 


for any hxed /c G M. This convergence, however, does not hold uniformly in k since Be is a singular 
perturbation of Bq. Using the uniform positivity of these symbol functions, we easily demonstrate 
the existence of the inverse operators 


B-\B^^ : L2(M)^L2(M), 

where Bq^ maps actually into the Sobolev space W^’^(M) and is hence smoothing because l/bo{k) 
decays quadratically at inhnity. The inverse of Be, however, is less regularizing because be{k) remains 
bounded as k ^ Too. In order to obtain asymptotic estimates for B~^, we introduce the cut-off 
operator 


Ue : L2(M) ^ L2(M) 


by defining its symbol function tt^ as follows 


TTeik) := 


1 for |A:| < - , 

0 else. 


(25) 


One of our key technical results is the following characterization of B~ ^, which reveals that Be admits 
an almost compact inverse. For m = 1, a similar but slightly stronger result has been given in [FP99, 
Corollary 3.5] using a careful Fourier-pole analysis of the involved integral operators. For m > 1, 
however, the symbol functions possess more poles in the complex plane and hence we argue differently. 




Figure 3: Sketch of the symbol function from (24), depicted on two intervals for e > 0 (black) and e = 0 
(gray). Notice that 6e(0) = inmi.^^bg{k) = 1 holds for all e > 0. 


Lemma 6 (asymptotic estimates for ^). For any e > 0, the operator respects the even-odd 
parity and is both self-adjoint and invertible on L^(M). Moreover, there exists a constant C such that 

\\U,B;^G\\^^^ + e-^\\{id-U,)Bf^G\\^<C\\G\\, (26) 

holds for all G G L^(M) and allO < s <1. Here, 11-112 2 denotes the usual norm in 


Proof In view of (22), (24) and Lemma 3, it remains to show (26). Using the properties of the sine 
function, see Figure 2, we readily demonstrate 


1 > 1 — sinc^ (mz) > 
Consequently, we get 


(min{|z| , 2})^ 
6 


for all z G M and m G N. 


1 — sinc"^ {mekf2) > — 


1 f for \k\ < f 


24 16 


else 


for all m, and hence 


be{k) > c 


1 + for |A;| < | 

1 /e^ else 


for some positive constant c > 0. Moreover, noting that 

G{k) 

be{k) 


B7^G{k) = 


and using Parseval’s theorem we estimate 


IWg||L=/' + + 


jgwi j,., 1 f 

b,(kf - cU|/.|<l l + 2fc" + 




as well as 




\k\>l beiky 


so (26) follows immediately. 


□ 


There exists another useful characterization of which relies on the non-expansive estimate 
PmelUlloo < ll^lloo> See Lemma 3. 
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Lemma 7 (von Neumann representation). We have 




a V 


i+l 


i=o (^e2 + J2m=i amm^) 

where the series on the right hand converges for any W £ (M). 

Proof. In the first step we regard all operators as defined on and taking values in L'^ 
use the abbreviation 


We also 


X, := 


Em=l 

£2 +eg 


and notice that (8) and (22) imply 


f 2 -£e2 


Since the operator norm of - computed with respect to the oo-norm ~ satishes 


T < _y_ 

ll^ellop-^2 + ,2 


< 1, 


the von Neumann formula provides 


B7^ = 


+ ci\ 


(id+X, 


-^e + 


£2 + ci 


id+- 


£2 + ci 


id +X^ + ... 


(27) 


in the sense of an absolutely convergent series of L°°-operators. In the second step we generalize this 
result using the estimates from Lemma 3. In particular, the right-hand side in (27) is well-defined 
for any W £ L2(M) since Lemma 3 ensures XgVL G L°°(M). □ 

Corollary 8 (invariance properties of Bf^). The operator Bf^ respects for both £ > 0 and £ = 0 the 
nonnegativity, the evenness, and the unimodality of functions. 


Proof. For £ > 0, all assertions follow from the representation formula in Lemma 7 and the cor¬ 
responding properties of the operators Ame, see Lemma 3. For £ = 0 we additionally employ the 
approximation results from Lemma 5 as well as the estimates from Lemma 6. □ 

Note that all results concerning Bf^ are intimately related to the supersonicity condition eg > eg. 
In a subsonic setting, one can still establish partial inversion formulas but the analysis is completely 
different, cf. [HMSZ13] for an application in a different context. 


3 Proof of the main result 

In view of the wave-speed scaling (8) and the fixed point formulation (13), the rescaled traveling 
wave problem consists in finding solutions £ L2(M) to the operator equation 

BeWe = , (28) 

where the linear operator Bg has been introduced in (22). Moreover, the nonlinear operators 

M M 

Qe[W] := l^mm^AmeiAmeWf , Ve[W] ■= ^ mAn,e^'^{me^An,eW) (29) 

m=l m=l 
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Figure 4: Potential energy {left panel) and phase diagram {right panel) for the nonlinear oscillator ODE 
(9) with coefficients (32), which determines the KdV wave Wq- There exists precisely one homoclinic orbit 
(solid black curve in the right panel) which corresponds to the solitary wave Wq. The closed loops inside the 
homoclinic orbits correspond to periodic KdV waves, see [FML14]. 


encode the quadratic and cubic nonlinearities, respectively, and are scaled such that the respective 
formal e-expansions involve nontrivial leading order terms. In particular, we have 




£—^0 
-^ 


Qo[l^] := 



(30) 


for any fixed W G L^(IR), see (20). Note also that (28) always admits the trivial solution We = 0. 

In what follows we solve the leading order problem to obtain the KdV wave Wq, transform (28) 
via the ansatz (10) into another fixed point equation, and employ the contraction mapping principle 
to prove the existence of a corrector I 4 G Lgyg^(M) for all sufficiently small e > 0. In [FP99], the last 
step has been solved using a operator-valued variant of the implicit function theorem. 


3.1 The leading order problem and the KdV wave 

Passing formally to limit e —)• 0 in (28), we obtain the leading order equation 

BoWo = Qo[Wo], (31) 


which is the ODE (9) with parameters 


di := 


12 

}^m=l OimlTT- 


d‘2 : = 


Em=l 


(32) 


In particular, the leading order problem is a planar Hamiltonian ODE with conserved quantity 
E = 2 (IP')^ + ^d 2 W^ — jdilK^ and admits precisely one homoclinic orbit as shown in Figure 4. 

Lemma 9 (linear and nonlinear leading-order problem). There exists a unique solution Wq G 
Lgygjj(M) to (31), which is moreover smooth, pointwise positive, and exponentially decaying. Moreover, 
the -kernel of the linear operator Cq with 


CoV := BoV - MoV, 


MoV := 2 


^ M 

E 

\m=l 




m 


\WoV 


(33) 


is simple and spanned by the odd function Wq. 

Proof. The existence and uniqueness of Hq follow from standard ODE arguments and the identity 
CqWq = 0 holds by construction. Moreover, the simplicity of the L^-kernel of the differential operator 
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£o can be proven by the following Wronski-type argument: Suppose for contradiction that Vi, V 2 £ 
L^(M) are two linearly independent kernel functions of Cq such that ti;( 0 ) 7 ^ 0 , where 


uj{x) := det 


Viix) 1/2 (x)\ 

Vlix) V'(x)J ■ 


The ODE C^Vi = 0 combined with Vi G L^(M) implies that Vi and Vl are continuous functions with 


\V{x)\ + \Vl{x)\ 


|x|^oo 


> 0 , 


and we conclude that ti;(x) —)• 0 as |a:| —)• 00 . On the other hand, we easily compute oj'{x) = 0 and 
obtain the desired contradiction. □ 

Since Wq is smooth, it satisfies (28) up to small error terms. In particular, the corresponding 
linear and the quadratic terms almost cancel due to (31). 

Lemma 10 (e-residual of Wo). There exists a constant C such that 

Qe[Wo]-B,Wo 


\Re\\2+\\Se\\2<C with Re := Se := Ve[Wo] 

holds for all 0 < e < 1. 

Proof. We first notice that Lemma 3 ensures 

\\Al,Wo\\, < \\A^eWo \\2 < IIW0II2 , WAmeWoW^ < ||Wo|L 

and in view of Assumption 1 we find 

M 


1 

ll'S'ell2 - \\ms‘^AmeWo\\l^ 11ILo11 2 < c. 


m=l 


Thanks to the smoothness of Wo, Lemma 5 provides a constant C such that 


AmsW^ - Wo' 


holds for j G {1,2}, and this implies 


+ 


AmeW^ - Wo' 




Ame(^AmeA^0^ Wq < AmeiyAmeW^Q^ Ams^^O 


+ 


< 


[AmeWof - 


+ Cm £ 


2 

2^2 


.A^^Wo" - Wo^ 


< (P™,Wo|L + ||Wo|L) P„^,Wo-Wo ||2 + CmV 

< Cm?£^ 


and hence 

||Qe[Wo]-Qo[Wo]||2 = 

Therefore, and since Wo satisfies (31), we get 

||S,Wo-SoWo|| 


M / M \ 

fdrnm^AmeiAmeWof - i ^ j Wq^ 

m=l \m=l J 


< Ce^ 


M 


We|l2 — 


+ c<E 




ALWo - Wo - ^W" 


12 '^'^0 


+ c, 


(34) 


(35) 


m=l 


12 



where the second inequality stems from the definitions of Be and Bq, see (22) and (23). Lemma 5 
also yields 


AmeWo - ILo - -^W!' 


24 


< Cm^e^ , 


AmeW” - W”\\^ < 


and combining this with (15)2 and the identity 


2 2 

Ai,Wo - hLo - = (A^e + id) (AmeWo - ILo - 


2^2 


me 


W” + 


24 


AmeWi;-Wi;\ , 


we arrive at 


AleWo -Wo- 


e'^m? 


24 


-VLn 


< 2 


2^2 


AmeWo - VLo - 


me 


-ILn 


+ 




^\\AmeW;,^-W(;\^ 


< Cm*e* 


The desired estimate for Re is now a direct consequence of (35). □ 

For completeness we mention that 

Wo{x) = ^ sech^(^dy^a;^ 

can be verified by direct calculations and that formulas for the spectrum of Cq can, for instance, be 
found in [MF53, page 768]; see also [FP99, Lemma 4.2]. 


3.2 The linearized traveling wave eqnation for e > 0 

For any e > 0, we define the linear operator Ce on L^(M) by 

M 

CeV := BeV - MeV , MeV := 2 Pmrn^Ame {(AmeWo) (AmeV)), (36) 

m=l 

where IFo G Lgygn(M) is the unique even KdV wave provided by Lemma 9. This operator appears 
naturally in the linearization (28) around Wo as 

Be{Wo + e^V)-Qe[Wo + e‘^V] = -e^Re + e^CeV - e^QeiV) 

holds due to the linearity of Be and the quadraticity of Qe- 

Lemma 11 (elementary properties of Ce). For any e > 0, the operator Ce is self-adjoint in L^(M) 
and respects the even-odd parity. Moreover, we have 

CeW c^w strongly in L2(M) 


for any W G W^’^(M). 

Proof Since IFo is smooth and even, all assertions follow directly from the properties of Ame and 
Be, see (22) and Lemma 5. □ 

Our perturbative approach requires to invert the operator Ce on the space Lg.^,gj^(M) - see the 
fixed point problem in Theorem 13 below ~ and in view of Lemma 9 one easily shows that Co has 
this properties. The singularly perturbed case e > 0, however, is more involved and addressed in 
the following theorem, which is actually the key asymptotic result in our approach. Notice that 
the analogue for M = 1 is not stated explicitly in [FP99] although it could be derived from the 
asymptotic estimates therein. 
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Theorem 12 (uniform invertibility of There exists 0 < e* < 1 sueh that for any 0 < e < e* the 
operator is continuously invertible on Lgyg^(M). More precisely, there exists a constant C which 
depends on e* but not on e such that 

\\Cf^G\\^<C\\G\\, 

holds for all D < e < and any G G Lgygj^(M). 

Proof. Preliminaries : Our strategy is to show the existence of a constant c* > 0 such that 

||/:eb^|l2>C* 11^112 (37) 

holds for all V G Lgygj^(M) and all sufficiently small e > 0, because this implies the desired result. In 
fact, (37) ensures that the operator 

: Levgn(M) —^ Lgygn(]R) 

has both trivial kernel and closed image. The symmetry of gives 

ker Ce = coker 

and due to the closed image we conclude that is not only injective but also surjective. Moreover, 
the e-uniform continuity of the inverse is a further consequence of (37). 

Now suppose for contradiction that such a constant c* does not exist. Then we can choose a 
sequence (en)„gN C (0, 1] with e^ 0 as well as sequences (I4)„gN ^ *-g^gn(M) and (Gn)„gN C 
Lgygjj(M) such that 

Ce„Vn = G^, ||K||2 = 1, l|Gn|l2 0. (38) 

Weak convergence to 0: By weak compactness we can assume that there exists Voo G Lgyg^(M) 
such that 


Vn Voo weakly in (M), 


(39) 


and using Lemma 11 we find 

(14o, dlocf)) = lim {Vn, = lim f) = lim (G„, <(>) = 0 

n^oQ n^oo n^oo 

for any sufficiently smooth test function f. In view of the definition of the differential operator Co - 
see (23) and (33) - we estimate 


Wo{x)(f)''{x)<lx <C'||i;i)||2 

R 

for all (j) G W^’^(M) and conclude that Voo belongs to W^’^(M), where CqVoo = 0 holds due to 

{CoVoo, 4>) = {Voo, c,o4>) = 0 • 

In other words, the even function Voo belongs to the kernel of Co and 


= 0 


follows from Lemma 9. 

Further notations: For the remaining considerations we abbreviate the constant from Lemma 6 
by D and denote by G any generic constant (whose value may change from line to line) that is 
independent of n and D. We further choose K > M sufficiently large such that 


sup ITo(0 < 
\i\>K-M 




, 771 '" 


( 40 ) 
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and denote by xk the characteristic function of the interval Ik ■= [—K, +K]- In what follows we 
write Vn = Vn^'^ + with 

Vhi) := XK , vP := (1 - xk) n.„K , := (id - U,jVn 


and observe that these definitions imply 


^ 11^^112 = 1 - 


* 6 { 1 , 2 , 3 } 


We also set 


i/h) ;= yh) 


(41) 


(42) 


and combine Lemma 3 with the smoothness of ILq to obtain 

l|[/h)|| <(7111/^)11 

II II 2 — ^ I wII 2 ' 

Moreover, by construction we have 

K = (uW + C/(2) + up + Gn) , 

SO the estimate 

+ yP\\2,2+^n^\H"\ ^ + I|f^^'^ll2 + l|f^F|l2 + l|G'n|l2) (43) 

is provided by Lemma 6. 

Strong convergence ofVp'^ and vjp: Inserting (38), (41), and (42) into (43) gives 


\\yp + yp\, + ep\\vp\ < CD (44) 

and hence 

vP 0 strongly in (M). (45) 

Thanks to 

vP=t) in L\Ik) (46) 

we also infer from (44) the estimate 

IkLIIw,, <IIW + rW||^^,<cD, 

where ||•||2 2 /^ denotes the norm in W^’^(/ii-). Since W^’^(/x) is compactly embedded into 1^(1;^), we 

conclude that the sequence (I4i^^)^gp^ is precompact in 1?{Ik)- On other hand, the weak convergence 
(39) combined with (45) and (46) implies 




Voo = 0 weakly in L^{Ik), 


and in summary we find vP —)■ 0 strongly in \J‘{Ik) by standard arguments. This even implies 

Y^) 0 strongly in L^(M) (47) 


as vanishes outside the interval Ik- 

(2^\ ^2') 

Upper bounds for \\Un q| 2 : Since the functions Vn are supported in R\ Ik, the functions 
f2) 

Ame„yn are supported in M \ lK-me„l 2 = {aJ : \x\ > K — men/‘i]. Moreover, we have 

|(^m£„ILo)(0| < sup ITo(x) < sup ITo(x) 

\x—^\<m£n/‘i \x—i\<M/2 
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for any given ^ G M. Therefore, and using 



for all X G M, 


we estimate 


VFo) m£n vP) 


< 


< 


sup |(“ 4 m£„lTo)(^)| j 

Ai\>K-M/2 I 


sup 

Ai\>K-M 




SO Lemma 3 gives 


and hence 


Ame„ {Ame„ Wq ) {Amen ^ ^ ^ 


( sup 

m>K-M 




( sup lTo(0 



1 

TD 


(48) 


due to (40) and (41). 

Derivation of the contradiction: Combining (43) with (42) gives 



and passing to the limit n —)■ oo we get 

limsup III4II2 ^ Ulimsup 11^41^^112 — \ 

n^oo n^oo 

thanks to (38)3, (45), (47), and (48). This, however, contradicts the normalization condition (88)2. 
In particular, we have shown the existence of a constant c* as in (37) and the proof is complete. □ 


3.3 Nonlinear fixed point argument 

Setting We = Wo + the nonlocal traveling wave equation (28) is equivalent to 

CeVe = Re+ Se + e^ Qe[K] + e" A4[K] 


with 

v.[r|:=Ll!L±fTC^. (4,, 

where Qe, Ve and Re, Se have been introduced in (29) and (34), respectively. Since Ce can be 
inverted for all sufficiently small e > 0, we finally arrive at the following result. 

Theorem 13 (existence and uniqueness of the corrector I4). There exist constants D > t) and 
0 < e* < 1 such that the nonlinear operator Te with 

Fe[V] :=Cf^(Re + Se + e^Qe[V]+e‘^Me[V]) (50) 

admits for any 0 < e < e* a unique fixed point Ve in the set Bd = {V : Lgygj^(M) : < D}. 
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Proof. Our strategy is to demonstrate that the operator Tg maps Bu contractively into itself provided 
that D is sufficiently large and e sufficiently small; the desired result is then a direct consequence 
of the Banach fixed-point theorem. Within this proof we denote by C any generic constant that is 
independent of D and e. We also observe that \Ar^Z\ < Ari iZj holds for any Z € and rj > 0, 

and recall that 


Re + SsH2<C 


is provided by Lemma 10. 

Estimates for the quadratic terms: For V £ Bd we find 

M / M \ 

|e'Qe[V^]| < e' E \\AmeV\\^Al,\V\ < , 

m=l \m=l / 

where we used the estimate (15)i, and in view of (15)2 we obtain 

|| e ^ Qe [ V ^]||2 < e ^ l’^CD ^ || y ||2 < e ^ l '^ CPA . 

In the same way we verify the estimate 

II M 

\e^Qe[V2]-e‘^Qem\\2< 


Y PmmH\\A^eV2\\^ + Pn^.Vl|loo ML|^2 " Vi | 


m=l 


< e^/^CD\\Al,\V 2 - Fillip < ||F 2 - V1II2 

for arbitrary Vi, V 2 £ B^. 

Estimates for the higher order terms: For Vi, V 2 G Bjo we set Zm,e,i '■= £‘^mAme{Wo + e'^Vi) and 
employ (15) 1 to estimate 


m,£,i||oo < WAmeWoW^ + \\AmeVi\ 
< £^m II Wolloo + ||y.||^ 


IZ, 

11 CXJ — 

< 

< £^m{c -h =: Cm,e ■ 

Due to the intermediate value theorem as well as the properties of we get 

M 


e^A 4 [h 2 ] - e^A4[Vi] < ^ m 




£ 


6 


m=l 

m=l 

^ m 'ymCm,e l“^m£F2 -4meLi| 
- 2^ ;^2 


< £ 


m=l 
2 


(c + £3/2^^ " ( 1 1/2 - Dl I 

\m=l / 


and hence 

\\e^Me[V2]-e^Me[Vi]\\^ < e2C7(c' + e3/22;)'||F2 - Vi|| 
after integration. A particular consequence is the estimate 

||A4[D]||2 < e‘^CD[c + 
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for any V £ Bd, where we used that A4[0] = 0. 

Concluding arguments: Combining all estimates derived so far with the dehnition of and the 
bounds for C~^ - see Lemma 12 - we verify 

||-^£[’C]||2 <C + e^/‘^CD‘^ + e‘^CD[c + 

for aWV £ Bo as well as 

\\Te[V2]\\,-\\Bem\,<(^e^/^CD + s‘^c[c + S^/^ l|V^2-Vi|l2 

for all Vi,V 2 G Bo- To complete the proof we first set D := 2C and choose afterwards e > 0 
sufficiently small. □ 

Corollary 14 (main result from §1). For any sufficiently small s > 0, the reformulated traveling 
wave equation (13) admits a unique even solution with speed such that 

\\We-Wo\\, + \\We-Wo\\^<Cs‘^ 

holds for some constant C independent of e. Moreover, We is nonnegative and smooth. 

Proof. The existence and local uniqueness of We = Wo + s'^Ve along with the L^-estimate is a direct 
consequence of Theorem 13. Moreover, re-inspecting the arguments from the proof of Theorem 13 
and using Lemma 7 we easily derive an uniform L°°-bound for the corrector Ve- Finally, the right 
hand side in (28) is - at least for sufficiently small e > 0 - nonnegative due to the properties of the 
KdV wave Wq and the potential <h, see Lemma 9 and Assumption 1. The nonnegativity of We is 
hence granted by Corollary 8. □ 

The constants in the proof of Theorem 13 are, of course, far from being optimal. In general, a 
solution branch s -£ We £ Lgyg^(M) on an interval [0, e*] can be continued for e > e* as long as 
the linearization of the traveling wave equation around We„ provides an operator Ce„ that can be 
inverted on the space Lgyg^(M). Since the shift symmetry always implies that W^^ is an odd kernel 
function of Ce„, the unique continuation can hence only fail if the eigenvalue of the linearized 
traveling wave operator 


M 

V^Yl m^A,^e.,K{melA^e.We^)Ame.,V (51) 

m=l 

is not simple anymore. Unfortunately, almost nothing is known about the spectral properties of the 
operator (51) for moderate values e*. It remains a challenging task to close this gap, especially since 
any result in this direction should have implications concerning the orbital stability of We„. 

For M = 1 it has also been shown in [FP99, Propositions 5.5 and 7.1] that the distance profile 
AeWe is unimodal (‘monotonic falloff’) and decays exponentially for x -£■ Too. For M > 1, it should 
be possible to apply a similar analysis to the velocity profile We but the technical details are much 
more involved. It remains open to identify alternative and more robust proof strategies. For instance, 
if one could show that the waves from Corollary 14 can be constructed by some variant of the abstract 
iteration scheme 


W^I3f^{Qe[W]+e^re[W]), 

the unimodality of We would be implied by the invariance properties of Ame and see Lemma 3 
and Corollary 8. A similar argument could be used for the exponential decay because Ame maps a 
function with decay rate A to a function that decays with rate 


A 


sinh (|emA) 
^emX 


18 



and since the von Neumann formula from Lemma 7 provides corresponding expressions for 
see [HRIO] for a similar argument to identify the decay rates of front-like traveling waves. In this 
context we further emphasize that only supersonic waves can be expected to decay exponentially. 
For subsonic waves with speed < Cq, the linearization of the traveling wave equation (3) predicts 
tails oscillations and hence non-decaying waves, see [HMSZ13] for a similar analysis with non-convex 
interaction potentials. 


List of symbols 


Pm 

linear and quadratic coefficients in force terms 

(5) 

7 m 

bounds for higher order force terms 

(6) 

Ce 

speed of the wave 

(8) 

Ue 

position profile 

(2) 

W, 

velocity profile 

(7) 

Ve 

corrector to the velocity profile 

(10) 

Re, Ss 

residual terms with respect to Wq 

(34) 

CLr; 

convolution operator and its symbol function 

(12), (16) 

Be, be 

auxiliary operator and its symbol function 

(22), (24) 

hie) ^e 

cut-off in Fourier space 

(25) 

Qe 

quadratic terms in We 

(29) 

Ve 

cubic and higher order terms in IF/ 

(29) 

Me 

quadratic combination of Wq and V/ 

(36) 

Be 

linear terms in V/ 

(36) 

Me 

remainder terms in V/ 

(49) 

Ve 

fixed point operator for Ve 

(50) 

Wq, Co 

velocity prohle and speed of the KdV wave 

Lemma 9 

Bo, Qo 

formal limits of Be, Qe 

(23), (30) 

Bo, Mo 

formal limits of Be, Me 

(33) 
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